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The Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope is a pair-conversion tele- 
scope designed to survey the gamma-ray sky from 20 MeV to several hundreds of GeV. In this energy 
band there are no astronomical sources with sufficiently well known and sharp spectral features to allow 
an absolute calibration of the LAT energy scale. However, the geomagnetic cutoff in the cosmic ray elec- 
tron-plus-positron (CRE) spectrum in low Earth orbit does provide such a spectral feature. The energy and 
spectral shape of this cutoff can be calculated with the aid of a numerical code tracing charged particles in 
the Earth's magnetic field. By comparing the cutoff value with that measured by the LAT in different geo- 
magnetic positions, we have obtained several calibration points between ~6 and ~13 GeV with an esti- 
mated uncertainty of ~2%. An energy calibration with such high accuracy reduces the systematic 
uncertainty in LAT measurements of, for example, the spectral cutoff in the emission from gamma ray 
pulsars. 

© 2011 Elsevier B.V. All rights reserved. 


1. Introduction 

The Large Area Telescope (LAT) is the primary instrument on 
the Fermi Gamma-ray Space Telescope, which was launched on 
June 11, 2008. The LAT is a 4 x 4 modular array of identical towers, 
each comprised of a tracker/converter (TKR) and a calorimeter 
(CAL), surrounded by an anti-coincidence detector (ACD). Each 
TKR module contains 18 layers of tungsten foil and paired silicon 
strip detector x-y tracking planes. Each CAL module contains a 
hodoscopic array of 96 Csl (Tl) scintillating crystals 8.6 radiation 
lengths deep. High energy gamma rays predominantly convert in 
the TKR into electron-positron pairs, and the energy of the devel- 
oping electromagnetic shower is measured in the CAL To cover the 
energy range of the instrument the individual CAL crystals must be 
able to accurately measure energy deposits between 2 MeV and 60 
GeV. To accomplish this, each CAL crystal can measure the depos- 
ited energy in four different overlapping energy ranges, namely 
2-100 MeV, 2 MeV-1 GeV, 30 MeV-7 GeV, and 30 MeV-70 GeV. 

The process of converting the scintillation light measured in 
each individual CAL crystal to an energy of the incident photon 
can be roughly divided into two phases. The first phase is the cal- 
culation of the energy deposition in each crystal (hereafter, the 
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“crystal energy”), which relies on calibration constants derived 
from the measured signals from sea-level muons and on-orbit cos- 
mic ray protons together with constants derived from an on-orbit 
electronic charge-injection system. The second phase is event 
reconstruction, which is comprised of a set of algorithms that ac- 
count for energy deposited in the TKR, leakage of the shower out 
the back and sides of the CAL, and energy lost in gaps between 
the towers. We developed the energy reconstruction algorithms 
[1] from a detailed instrument model and Monte Carlo simulation 
based on GEANT4. 

To verify the fidelity of the instrument model and Monte Carlo 
simulation, we performed an extensive beam test campaign in 
2006 at CERN and the GSI heavy ion accelerator laboratory. To 
eliminate the handling and shipping risks to the flight detector 
subsystems, and to reduce cost and schedule impacts, these tests 
were not done on the full LAT but on a Calibration Unit (CU), which 
consisted of two fully-populated flight-spare TKR modules and 
three fully-populated flight-spare CAL modules. A more detailed 
description of the beam test as well as the energy reconstruction 
algorithms will be presented in Section 2. 

To perform an in-flight verification of the absolute energy scale 
of the LAT it is necessary to find an astrophysical source with a 
spectral feature whose absolute energy peak and shape are well 
known. Because sources in the LAT energy band exhibit smooth 
spectral forms, it is not trivial to satisfy this requirement. Further- 
more, the LAT is far more sensitive than any previous gamma-ray 
telescope, so it is not possible to use existing measurements of 
celestial sources to perform this calibration. The LAT has shown 
that gamma-ray pulsars typically exhibit spectral cutoffs in the 
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1-10 GeV range [2] and, for bright pulsars, the spectral cutoff en- 
ergy can be measured with statistical uncertainty approaching 1% 
[3], If the systematic uncertainty on the energy scale in this energy 
range can be reduced to the few-percent level, precise measure- 
ments of spectral cutoff can address detailed questions of the emis- 
sion physics and geometry in pulsars [3], In the orbit of the Fermi 
observatory (inclination of 25.6° and altitude of 565 km), the geo- 
magnetic cutoff in the cosmic-ray electron plus positron (CRE) 
spectrum can serve as such a source in the range between ~6 
and ~13 GeV. Electrons and positrons are a good calibration source 
because they interact in the LAT CAL in the same way that photons 
do: they both produce electromagnetic showers. The energy mea- 
surement bias coming from the slight differences between CRE and 
photons interactions in the LAT is much less than 1% [4]. The capa- 
bilities of the LAT to measure the CRE spectrum have been well 
demonstrated [4,5]. The cutoff rigidity 4 can be predicted by numer- 
ically tracing particle trajectories in the Earth’s magnetic field, and 
the comparison of the predicted and measured values provides the 
opportunity to perform this validation. 

In order to measure the deviation between the reconstructed 
electron spectrum and the calculated value as a function of energy, 
this analysis is performed in several Mcllwain L intervals. The Mcll- 
wain L-parameter [6] is a parameter describing a set of the Earth’s 
magnetic field lines, in particular those which cross the Earth’s 
magnetic equator at a number of Earth-radii equal to the L-value. 
Magnetically equivalent positions (from the standpoint of the 
incoming charged particle) around the world will by definition 
have the same Mcllwain L values, therefore making this parameter 
particularly convenient for characterizing cutoff rigidities [7]. The 
orbital inclination of the Fermi orbit fixes the range of Mcllwain 
L values accessible for this analysis to 1.00-1.72. 


2. LAT energy calibration 

The absolute energy scale of the LAT is defined by comparing 
the signals in the CAL crystals with the amount of energy a Monte 
Carlo simulation indicates should be deposited by on-orbit relativ- 
istic protons. The details of the energy calibration, referred to as 
the "proton inter-range calibration”, are given in [8], To perform 
this calibration, we first correct the observed CAL signals for elec- 
tronic non-linearities (as measured by an electronic charge-injec- 
tion process) and for position-dependent scintillation response 
(as measured by a direct calibration with sea-level muons and 
on-orbit with protons). We then compare the corrected CAL signals 
with the distribution of deposited energies predicted by a GEANT4 
simulation of the on-orbit spectrum of primary cosmic-ray protons 
passing through the LAT. This calibrates the highest-gain, lowest- 
energy range of the CAL readout. The remaining gain ranges [1] 
of the CAL readout are calibrated by enforcing that adjacent gain 
ranges give the same measured energy in the regions of energy 
space in which they overlap [8], 

The Galactic cosmic-ray (GCR) element abundance peaks could 
in principle be used to calibrate the higher energy ranges of the 
LAT; however, the scintillation efficiency of heavy ions in Csl (Tl) 
is not identical to that of electromagnetic showers by an ele- 
ment-dependent factor that is not well known or measured. We 
measured the relative scintillation efficiencies of sub-relativistic 
ions in Csl (Tl) crystals in a series of beam tests [9], but we do 
not have sufficient confidence that these efficiencies apply at rela- 
tivistic energies. Thus, we are unable to relate the observed signals 


4 Rigidity is defined as the particle momentum divided by its charge. The cutoff we 
are measuring in the CRE spectrum is not the vertical cutoff but rather averaged over 
all angles. 


from GCR element peaks to electromagnetic shower energy depo- 
sitions accurately enough to calibrate the LAT energy scale. 

The geomagnetic cutoff in the CRE spectrum (measured in the 
energy range ~6 to ~13 GeV) is a good alternative source to cali- 
brate the higher energy ranges of the CAL readout. At these ener- 
gies, the maximum energy per crystal is of the order of ~1 GeV, 
and will be read out in the higher energy ranges of the CAL. As a 
reference, 1 GeV of energy per crystal is equivalent to 100 times 
what a minimum ionizing particle releases per crystal. 

The LAT CU was also calibrated using sea-level cosmic-ray 
muons, but the absolute energy scale was then cross-checked 
against calibrated beam lines at CERN, showing an energy deposit 
systematically larger than expected. Such a direct calibration relies 
on the knowledge of the incoming beam energy (determined with 
a 1% accuracy), the geometry of the beam line and the CU detector, 
the Monte Carlo of the electromagnetic shower development with- 
in the detector, as well as on a good control of the environmental 
effects that might influence the CU response (temperature, humid- 
ity, exceedingly high particle rates). 

The energy reconstruction of the LAT (as well as of the CU) is 
based on three different algorithms: a parametric correction based 
on the barycenter of the shower, a fit to the shower profile taking 
into account the longitudinal and transverse development of the 
shower, and a maximum likelihood fit based on the correlation be- 
tween the total deposited energy, the energy deposited in the last 
layer of the CAL and the number of tracker hits. For each event, the 
best energy reconstruction method is selected by means of a clas- 
sification tree analysis described in [1], We would like to empha- 
size that at the energies considered in this analysis (i.e. between 
~6 and ~13 GeV) the correction factors for losses due to leakage 
out of the detector are reasonably small (on average of the order 
of 30%) when compared to higher energies (of the order of 50% 
for 100’s of GeV) and therefore the method described here is test- 
ing a mixture of the crystal energy calibration and the leakage 
corrections. 

Following the beam test campaign, the GEANT4 simulation of 
the CU was updated to reflect the best available description of 
the CU detector, the beam line, and the particle interactions. Most 
notably, the routines describing the Landau-Pomeranchuk-Migdal 
effect (LPM [10]) in electromagnetic shower developments were 
updated after finding that initial GEANT4 implementations were 
not providing a satisfactory description of the shower longitudinal 
development. The energy resolution for all three available energy 
reconstruction methods was measured to be consistent with 
expectations. However despite these improvements, the recon- 
structed energy in the CU was consistently higher than the beam 
energy by ~9%, on average, with further fluctuations of ~5%, 
depending on incoming beam energy (5-282 GeV), angle (0-60 de- 
grees) and position of the crystal within the shower [4,11]. 

The energy scale factor derived from the direct CU calibration 
was not applied to the LAT because of the differences in the beam 
test and on-orbit environment, particularly the distribution and 
rate of the particles incident on the calorimeter, and the tempera- 
ture variations for the CU during the data taking, which were con- 
sidered responsible for the different LAT and CU energy scales. 


3. Particle tracing 

Geomagnetic cutoff rigidities can be obtained by tracing cos- 
mic-ray trajectories in a model of the Earth’s magnetic field. The 
standard mathematical description of the Earth's magnetic field 
is given by the International Geomagnetic Reference Field models 
(IGRF) [12], These models consist of the Gauss coefficients defining 
the spherical harmonic expansion of the geomagnetic potential up 
to a given order. In the IGRF-1 1 model, used for this analysis, the 
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Fig. 1. Illustration of e and e* of various energies traced in the Earth’s magnetic field from a point at longitude and latitude (0,0) and 565 km above the Earth’s surface 
looking down on the North Pole. The vertical cutoff is taken to be the lowest rigidity which the particle needs to be allowed when traced vertically through the Earth’s 
magnetic field. As the rigidity of the particle decreases the amount of geomagnetic bending increases and when the trajectory intersects the Earth it is taken to be forbidden 
(i.e. of secondary origin). The cutoff values are different for positively and negatively charged particles because the Earth's magnetic field is not a perfect dipole. Particles 
labeled as allowed are taken to be of Galactic origin. Earth and trajectories are to scale. 


maximum multipole moment of the expansion is the 13th. The 
IGRF model is computed by the participating members of the inter- 
national association of geomagnetism and aeronomy (IAGA) work- 
ing group V-MOD [12] which is supported by the organizations 
involved in operating magnetic survey satellites, observatories, 
magnetic survey programs and World Data Centers. IGRF-1 1 pro- 
vides a definitive main field model for epoch 2005, a main field 
model for 2010, and a linear predictive secular variation model 
for 2010-2015. We use the particle trajectory tracing code (hereaf- 
ter tracer) developed by Smart and Shea [7] and the 1GRF-11 in this 
analysis to measure the geomagnetic cutoff. 

Since it is difficult to trace the trajectory of an incoming particle 
through the magnetic field and expect it to intersect the precise 
location desired, it is more efficient to calculate the trajectory in 
the reverse direction (i.e. trace out a positron to mimic an incoming 
electron). Therefore the starting point of the trajectory is given by 
the geographic coordinates, altitude and orientation of the space- 
craft and the trajectories are propagated using IGRF-1 1 [12]. 

At each of a grid of locations covered by the Fermi orbit, evalu- 
ated every 15 s from Aug 2008 through Jul 2009, we simulated an 
ensemble of test particles. We selected energies of the test parti- 
cles according to the power-law spectrum measured by the LAT 
[4], and generated both electrons and positrons in the abundance 
ratio measured by PAMELA [13]. We use the tracer code to deter- 
mine for each test particle whether it could have originated from 
outside the geomagnetosphere. Trajectories that eventually inter- 
sect the Earth’s atmosphere 5 are rejected and labeled as forbidden 
(secondaries hereafter). Trajectories that reach 20 Earth radii are ac- 
cepted as allowed (or of Galactic origin, primaries hereafter). Fig. 1 
illustrates a few sample trajectories, both forbidden and allowed. 

We selected electron and positron data from the first year of 
Fermi LAT data taking (Aug 2008 through Jul 2009) using the meth- 
ods described in [4] for the 100 MeV to 100 GeV energy range. In 
Fig. 2 are the distributions (for E > 30 GeV) of the reconstructed an- 
gle with respect to local zenith (upper panel) and azimuth (lower 
panel) for the flight data and for the tracer output. The overall 


5 Distance taken to be 20 km from the Earth’s surface. 




Azimuth 

Fig. 2. The reconstructed angle with respect to local zenith (upper panel) and 
azimuth (lower panel) for data and tracer. In the lower panel. North is at 0°, East at 
90°, South at 1 80° and West at 270°. The LAT acceptance for electrons and positrons 
as well as the rocking profile have been convolved in these distributions and 
influence their shape. Both are averaged over the orbit and for energies greater than 
30 GeV in order to compare the distribution of primary CREs. There is an overall 
good agreement between data and tracer. 


agreement between flight data and the tracer output is very good, 
demonstrating that the angular distributions have been well 
described. 
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4. Analysis 

4.J. Estimating the fraction of primary cosmic ray electrons 

The population of CREs in low Earth orbit is a mixture of pri- 
mary and secondary cosmic rays, where primary implies of extra- 
terrestrial origin and secondary implies a mixture of splash and re- 
entrant particles produced in the interaction of primaries in the 
Earth’s atmosphere. The trajectories of the secondary particles 
are extremely difficult to simulate reliably. By definition our trac- 
ing code only provides information on the allowed (i.e. primary) 
particles. As a consequence it is necessary to estimate the fraction 
of secondary particles, as a function of energy, from the flight data. 
This constitutes one of the most delicate aspects of the analysis be- 
cause the cutoff rigidity is found by fitting the spectrum (as will be 
described in detail in Section 4.2) and its shape is influenced by the 
fraction of secondaries. 

The azimuthal distribution in Earth centered coordinates of the 
secondary population is different from that of the primary parti- 
cles. We can exploit this fact and perform a template fit to identify 
the fraction of each population. To perform this task it is first nec- 
essary to choose the appropriate templates to describe the popula- 
tions. It is safe to assume that the population at low energy in flight 
data (E < E c , where E c is the energy corresponding to the cutoff 
rigidity) is predominately composed of secondary particles, and 
can therefore be used as a template. The output from the tracer 
code provides the template for the primaries. Fig. 3 shows an 
example of the template fitting performed to estimate the fraction 
of primaries in the interval 1.0 < Mcllwain L<1.14. The gradual 
transition from a pure sample of secondary particles (top left panel, 
for energies between 2 and 4 GeV) to that of only primaries (bot- 
tom right panel, for energies between 16 and 18 GeV) is evident. 


§ 1 

o 

TO 

LI- 

S' 

TO 

E 

□I 0.5 


0 


Fig. 4. Fraction of primary electrons and positrons as a function of energy for the 
Mcllwain L interval 1.00 <L< 1.14. Note that the primary fraction does not reach 
100% because contamination from hadrons has not yet been accounted for. This is 
done when reconstructing the primary spectrum as described in Section 4.2. 


10 4 


Energy (MeV) 


For reference, the cutoff energy in this Mcllwain L interval is 
~13 GeV. For simplicity, only four energy intervals are shown. 
However the analysis has been performed in 16 overlapping en- 
ergy intervals. The resulting value for the fraction of primaries as 
a function of energy for this same Mcllwain L interval is shown 
in Fig. 4. 

4.2. Measuring the cutoff rigidity 

The primary e~ + e + spectrum can be parameterized by [4]: 

dN(E) cE- r 

dE (1 + (E/E c )- 6 ) 




0 50 100 150 200 250 


300 350 

Azimuth 




Fig. 3. Azimuthal distribution for flight data, tracer and secondaries in the Mcllwain L interval 1 .0 < L < 1 .14 averaged over all zenith angles. The black distribution is the linear 
combination of tracer and secondaries which best fits the data. Each panel depicts the distributions in a given energy interval (labeled in the panel). In the upper left panel, 
where the secondary population is dominant, a maximum is evident in the westward direction where the secondary positrons are dominant. With increasing energy, the 
primary population gradually increases in the eastward direction. The overall shape of the azimuth distribution is affected by the LAT acceptance and the rocking profile. The 
cutoff rigidity in this Mcllwain L range is ~13 GeV. The fraction of primaries evaluated from this template fitting is shown in Fig. 4. 
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where r is the spectral index and E c is the energy corresponding to 
the cutoff rigidity. We fit the primary CRE spectrum (both data and 
tracer) with Eq. (1) to get the value of the cutoff rigidity. 

To reconstruct the CRE primary spectrum it is first necessary to 
remove the secondary population from the count spectrum and 
then correct for the hadronic contamination. The residual contam- 
ination from hadrons, h, is estimated by applying the CRE event 
selection to the on-orbit simulation used by the Fermi LAT collab- 
oration. We found that this residual hadron contamination ranges 
between ~15% below the cutoff and ~5% above with an estimated 
uncertainty on the absolute hadron flux of 20%. The on-orbit sim- 
ulation is based on data from CR experiments and includes all 
the components of the Galactic cosmic rays as well as the re-en- 
trant and splash Earth albedo particles. A detailed description of 
this simulation can be found in [4]. The residual contamination 
due to secondaries, s, is estimated as described in Section 4.1. 
The count spectrum is multiplied by the purity factor defined as 
P- F= (1 — h) ■ (1 - s) in each energy bin. This background-sub- 
tracted count spectrum is divided by the effective geometry factor 
(EGF) [4] and the width of the energy interval, to obtain the final 
spectrum. 

The tracer count spectra are reconstructed by requiring the al- 
lowed condition and dividing by the width of the energy interval. 
To properly compare the cutoff rigidities, the tracer spectra have 
been binned in measured energy (i.e., they have been convolved 
using the energy resolution of the Fermi LAT for electrons and pos- 
itrons [4]). No efficiency correction is required for the tracer spec- 
tra. To measure the spectra in Mcllwain L intervals an extra cut 
specifying the interval is applied to both data and tracer. 

A comparison between the counts spectra measured from flight 
data and from the tracer output is shown in Fig. 5. The shapes of 
the spectra are in good agreement, indicated by the ratio of the 
two depicted in the lower panel of the figure. Both the reconstruc- 
tion and fitting procedures are applied to each Mcllwain L interval 
considered in this analysis. 

4.3. Assessment of the systematic uncertainties 

The main sources of systematic uncertainty in this analysis are: 
the choice of the energy interval in which to fit the spectrum, the 
choice of the energy interval for the secondary template, the 
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Fig. 5. Energy spectrum of data and tracer for the Mcllwain L interval 1.0 <L< 1.14. 
The black line is a fit to the data using Eq. (1 ) and the resulting cutoff energy has a 
value of 13.27 ± 0.10 GeV, and is indicated by the dot dashed line in the figure. The 
lower panel depicts the ratio of the two spectra. 


accuracy of the geomagnetic field model, and the accuracy of the 
Monte Carlo simulations used for the estimation of the residual 
hadron contamination as well as the evaluation of the EGF. 

To estimate the systematic uncertainty from the choice of en- 
ergy interval used for the spectral fits, we performed a series of fits 
over a range of energy intervals narrower and broader than the 
chosen interval and calculated the root-mean square (rms) of the 
distribution of resulting values of the cutoff energy. The rms was 
found to be no greater than 0.8% for all the Mcllwain L intervals. 

Since the key ingredients in the calculation of the fraction of 
CRE primaries are the secondary and primary templates, it is cru- 
cial to validate them. The tracer output (used as the template for 
the primaries) can be validated directly with the Fermi LAT flight 
data, as was described in Section 3. The secondary template is 
more difficult to validate given that the output from the tracer 
code does not provide any information regarding the secondaiy 
populations. As stated in Section 4.1, it is safe to assume that the 
population at low energy in flight data ( E < E c ) is predominantly 
composed of secondary particles. However the choice of the energy 
interval is not well determined, and can therefore be a source of 
uncertainty. This choice can in fact influence the shape of the spec- 
trum and thus the final value for the cutoff rigidity. To investigate 
its effect, we chose several different energies ranges in which to 
define our secondary templates (for each Mcllwain L interval) by 
varying both the width of the energy interval as well as the dis- 
tance from the cutoff energy. With the resulting fraction of prima- 
ries obtained from each template we reconstructed the spectrum 
and obtained the cutoff rigidity following the procedure described 
in Section 4.2. The spread of the ratio of the cutoffs (found to be no 
greater than ~2% for all Mcllwain L intervals) gives an estimate for 
the uncertainty due to the secondaiy template. 

The accuracy of the cutoffs derived from the particle tracing 
code is limited by the uncertainty of the geomagnetic field model. 
Several cross checks on the accuracy of the predicted rigidity cut- 
offs have been performed using satellite experiments, in particular 
the cosmic-ray isotope experiment HEAO-3 C2 [14] (C2) and the 
Mass Spectrometer Telescope (MAST) on the Solar Anomalous 
and Magnetic Particle Explorer (SAMPEX) spacecraft [15], The 
HEAO-3 satellite flew between 1979 and 1981 at an altitude of 
496 km and inclination of 43.6°. From the comparison between 
the experimental and computed cutoffs of oxygen nuclei in the 
5 GeV/n range it was found that the calculated cutoffs were sys- 
tematically ~3-5±2% higher than the measured ones [14], The 
SAMPEX mission was launched into an orbit with 82° inclination 
with an average altitude of ~600 km and operated from 1992 to 
mid 2004. They performed a comparison between measured cutoff 
rigidities of 0.3-1. 7 GeV protons during geomagnetically quiet 
times with those predicted by the Smart and Shea tracer code 
(based on IGRF-10) and found that the latter systematically overes- 
timate the data by 8-14% [16]. 

Directly comparing these findings with the measurements pre- 
sented here is difficult because in both the SAMPEX and HEAO-3 
cases the measured cutoff rigidities are much lower than those 
appropriate for Fermi. Smart and Shea in [17] assert that the accu- 
racy of the model improves in regions closer to the geomagnetic 
equator (i.e. for larger cutoff rigidities). According to their analysis, 
it is reasonable to assume that the bias in the tracer output used in 
this analysis is no greater than ~3-5% for the Mcllwain L regions 
spanned by the Fermi observatory. However, this is just an extrap- 
olation based on the measurements performed by HEAO-3 and 
SAMPEX at lower rigidities and a definitive answer on the value 
of this bias at the Fermi rigidities remains unknown. 

The uncertainty of the hadronic component of the Monte Carlo 
simulations used to evaluate the residual contamination is esti- 
mated to be no greater than 20%. Its contribution to the overall sys- 
tematic uncertainty in the final measurement of the cutoff is small 
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Fig. 6. The path length-corrected crystal energies of the B,C,N, and O peaks 
corresponding to the first four months of Fermi LAT operation. The large apparent 
abundance of B (compared to secondary and primary Galactic cosmic-ray 
abundance) is due to secondary B produced by primaiy C interacting in the ACD 
and TKR. 

(0.5%) because the contamination is always below 20% in every 
Mcllwain L interval. The assessment of the systematic uncertainty 
due to the EGF is described in detail in [4] and for this energy range 
has a value of 10-15%. We studied its effect on the final value of the 
cutoff via a bracketing method. To do this we generated a set of 
EGFs that have a maximal impact on the final shape of the spec- 
trum, reconstructed the spectrum using these new EGFs and ana- 
lyzed the final value of the cutoff energy obtained with these 
new EGFs. The resulting values varied by ±1%. 

4.4. Cross check of method 

To demonstrate that this analysis has the sensitivity to measure 
deviations of several percent in the energy scale of the LAT, we re- 
peated the analysis on a test sample of traced particles with an 
added +5% shift in energy. The resulting measurement of the cutoff 
energy was 5.4 ± 0.2% higher, which is consistent with the bias we 
introduced in the input. Thus we conclude that indeed this method 
is capable of measuring energy scale errors of the magnitude we 
expect. 

Pre-launch tests of crystal boule samples showed [18] that the 
light yield of the Csl (Tl) crystals in the LAT calorimeter could be 
expected to decrease by ~1% per year from radiation damage in 
the charged-particle environment of low Earth orbit, primarily 
from trapped particles in the South Atlantic Anomaly. 

By analyzing the path length-corrected 6 crystal energies of four 
abundant GCR elements, namely boron, carbon, nitrogen and oxy- 
gen, over time it is possible to verify this prediction with flight data. 
In fact, all four of these peaks show a similar linear decrease per year. 
These results are consistent with the pre-launch predictions. An 
example of these peaks is shown in Fig. 6. As a further cross check 
of the sensitivity of the method used in this analysis, we have com- 
pared the cutoff rigidity values measured in the first 60 days to the 
last 60 days of the first year of operations. We found that the mea- 
sured cutoff energy decreased by 1.9 ± 0.9% over this time interval, 
which is consistent with the decrease from radiation damage mea- 
sured with the GCR element peaks, as can be seen in Table 1. We 
note also that this provides further confirmation of the sensitivity 
of our method. 


6 The energy deposit is normalized by cosd, where 0 is the incidence angle of the 
particle. 


Table 1 

Drift per year for four abundant GCR element (B,C,N, and O) 
peaks and the CRE cutoff energy measured in Fermi LAT flight 
data over the first year of operations. 


Drift (%/yr) 

B 

1.49 ±0.05 

C 

1.60 ±0.04 

N 

1.50 ±0.06 

O 

1.46 ±0.03 

CRE 

1.90 ±0.90 



Fig. 7. Distribution of cutoff energy ratios as a function of Mcllwain L. The ratios are 
constant within the errors over the entire range. The quoted value for the ratio is 
the weighted mean. The statistical errors are represented by the black error bars, 
and the systematic errors are in gray. As a reference, the value for the cutoff energy 
in each Mcllwain L interval is also shown. 

5. Results 

We applied this analysis to a data set spanning the first year of 
LAT operations. We found that the measured cutoff energy exceeded 
the predicted cutoff energy by 1 .026 ± 0.005 (stat) ±0.025 (sys) in the 
6 GeV to 1 3 GeV range. The systematic error is the sum in quadrature 
of our estimates in Section 4.3, excluding those from the uncertainty 
in the magnetic field model and particle tracing code. As shown in 
Fig. 7, the ratio between measured cutoff and predicted cutoff is con- 
stant within the errors over this energy range. 

For the energy range in question— where the instrument accep- 
tance is relatively flat— an energy-independent bias in the absolute 
energy scale translates into a rigid shift of the spectrum, scaling 
the normalization but leaving the spectral index unaffected. For 
the physically motivated case of a source with a power-law spec- 
trum with index r = — 2, the shift in the normalization is of the same 
amount as the bias. Therefore the effect of an error at the ~2% level in 
the energy scale is smaller than other instrument-related sources of 
systematics (most notably the uncertainties in the effective area). 

6. Conclusions 

We provide a method for calibrating the LAT energy scale based 
on the measurement of the geomagnetic cutoff in the primary CRE 
spectrum. The results shown here confirm that the LAT calorimeter 
calibration procedure, based on measured electronic gains and 
non-linearities, observed proton signals, and enforced consistency 
between readout energy ranges, is accurate to the ~2% level. It is 
worth noting that the method developed here is the only such cal- 
ibration which relies uniquely on flight data and uses electromag- 
netic showers in the LAT. Accordingly, this is a very important cross 
check on other methods such as beam tests and calibrations with 
cosmic-ray heavy ions. 

Our conclusion relies, however, on the fidelity of the current IGRF 
and particle trajectory tracing code. If we take the measurements of 
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HEAO-3 C2 and SAMPEX to indicate a bias of the order of +3-5%, the 
absolute energy scale of the LAT may have to be decreased by 5-7%. 
This conclusion, if confirmed, would make the energy scale more 
consistent with the beam test results and the Monte Carlo simula- 
tions of the CU. This would suggest that the on-orbit energy calibra- 
tion technique for the LAT should be revised. However, the cross 
checks reported in [ 14,1 6] on the predictions of the tracer code were 
both performed where the cutoff rigidities are much lower than 
those measured by the LAT and in geomagnetic regions where the 
trajectory tracing calculations are less accurate [7 ]. It is therefore un- 
clear whether this bias in the theoretical cutoffs also applies to the 
results presented in this work. We therefore prefer to determine 
the energy scale of the LAT using the results of the proton inter-range 
calibration, which are confirmed by the independent measurements 
discussed in this paper. 
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